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Abstract 

We address the generation of fully inseparable three-mode entangled states of 
radiation by interlinked nonlinear interactions in x media. We show how 
three-mode entanglement can be used to realize symmetric and asymmetric 
telecloning machines, which achieve optimal fidelity for coherent states. An 
experimental implementation involving a single nonlinear crystal where the two 
interactions take place simultaneously is suggested. Preliminary experimental 
results showing the feasibility and the effectiveness of the interaction scheme 
with seeded crystal are also presented. 
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1. Introduction 

The successful demonstration of continuous variable (CV) quantum teleportation^* 3 - and 
dense coding 4 * opened new perspectives to quantum information technology based on Gaus- 
sian states of light. Besides having been recognized as the essential resource for teleportation 1 " 
and dense coding 5 -, the entanglement between two modes of light has been proved as a valu- 
able resource also for cryptography" improvement of optical resolution 8 , spectroscopy^, 
interferometryiS, state engineering^, and tomography of states and operation ^ 12 ! 13 . 

These achievements stimulated a novel interest in the generation and application of multi- 
partite entanglemen t 14 i 15 i 16 i 17 , which has already received attention in the domain of discrete 
variables. Multipartite CV entanglement has been proposed to realize cloning at distance 
(telecloning) 18 ! 1 ? , and to improve discrimination of quantum operations^. The separability 
properties of CV tripartite Gaussian states have been analyzed irA where they have been 
classified into five different classes according to positivity of the three partial transposes 
that can be constructed. Moreover, it has been pointed out that genuine applications of 
three-mode entanglement requires fully inseparable tripartite entangled states 2 ^, i.e. states 
that are inseparable with respect to any grouping of the modes. 

Experimental schemes to generate multimode entangled states have been already sug- 
gested and demonstrated. The first example, although no specific analysis was made on 
the entanglement properties (besides verification of teleportation), is provided by the orig- 
inal teleportation experiments of Ref.— where one party of a twin-beam (TWB) was mixed 
with a coherent state. A similar scheme, where one party of a TWB is mixed with the 
vacuum^ 5 - has been demonstrated, and applied to controlled dense coding. More recently, 
a fully inseparable three-mode entangled state has been generated and verified^ by mixing 
three independent squeezed vacuum states in a network of beam splitters. In addition, a 
four-mode entangled state to realize entanglement swapping with pulsed beams have been 
generated^ 3 -. 

All the above schemes are based on parametric sources, either of single-mode squeez- 
ing or of two-mode entanglement i.e. TWB, with multipartite entanglement resulting from 
further interactions in linear optical elements (e.g. beam splitters). In this paper, we fo- 
cus on a scheme involving a single nonlinear crystal, in which the three-mode entangled 
state is produced by two type I, non-collinearly phase-matched interlinked bilinear interac- 
tions that simultaneously couple the three modes 2 ^. A similar interaction scheme, though 
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realized in type II collinear phase-matching conditions, is described in Ref.— . Compared 
to this work, our choice of non-collinear phase-matching provides remarkable flexibility to 
our experimental setup, whereas the choice of type I interaction prevents the generation of 
additional parties. Moreover, we avoid the losses brought about by the mode-matching in 
multiple beam splitters in that we achieve the three-partite entanglement as soon as we find 
the configuration that fulfills the phase-matching condition for both interactions. 

The paper is structured as follows. In Section El we describe the generation of three-mode 
entanglement in a single nonlinear crystal where two interlinked bilinear interactions take 
place simultaneously. We obtain the explicit form in the Fock basis of the outgoing three- 
mode entangled state, and also address the characterization of entanglement. In Sections 
E] and E] we show how the three-mode entangled state obtained in our scheme, either for 
initial vacuum state or by seeding the crystal, can be used to build symmetric and asym- 
metric telecloning machines that achieve optimal fidelity for coherent states. In Section |3] 
we show how three-mode entanglement may be used for conditional generation of two-mode 
entanglement, in particular of TWB state. The scheme is of course less efficient than di- 
rect generation of TWB in a parametric amplifier, but it may be of interest in applications 
where entanglement on-demand is required. In Section 13 we discuss the experimental im- 
plementation of our generation scheme. We show the feasibility of experiments in the case 
of interaction with seeded crystal and report preliminary experimental results. Section [7| 
closes the paper with some concluding remarks. 

2. Generation of three-mode entanglement 

The interaction Hamiltonian we are going to consider is given by 



Hi n t describes two interlinked bilinear interactions taking place among three modes of the 
radiation field. It can be realized in x media by a suitable configuration which will be 
discussed in Section El The effective coupling constants jj, j = 1,2, of the two parametric 
processes are proportional to the nonlinear susceptibilities and the pump intensities. The 
Hamiltonian in Eq. (jl} has been firstly studied in2£ though not for the generation of 
entanglement. The Hamiltonian admits the following constant of motion 



Hi 



'int = 7l«l a 3 + 72<4 a 3 + h.C. . 



(1) 



A(i) = N x (t) - N 2 (t) - N 3 (t) = A(0) , 



(2) 
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where Nj(t) = (aAt)a(t)) represent the average number of photons in the j-th mode. If we 
take the vacuum |0) = |0)i <8> 1 0) 2 <£> 1 0)3 as the initial state we have A = i.e. Ni(t) = 
N 2 (t) + N 3 (t) Vi. The expressions for Nj(t) can be obtained by the Heisenberg evolution of 
the field operators, which read as follows 

a\(t) = /ial(0)+/ 2 a 2 (0) + /3a 3 (0) 
a 2 (t) = gia\(0) + g 2 a 2 (0) + g 3 a 3 (0) 

a 3 (t) = hxaliO) + h 2 a 2 (0) + h 3 a 3 (0) . (3) 
The explicit expressions of the coefficients fj, gj and hj, j = 1, 2, 3, are obtained in appendix 



1X1 we omit the time dependence for brevity. By introducing Q = a/|7 2 | 2 — |7i| 2 we have 

Ni = N 2 + N 3 , 

N 2 = J^&!![oosJtt-l] a , 

iV 3 = ^ sin 2 (fit). (4) 
The evolved state reads as follows 27 



(5) 



where U t = exp (—iH int t) is the evolution operator, and we have already used the conserva- 
tion law. The state in Eq. (JHJ) is Gaussian, as it can be easily demonstrated by evaluating 
the characteristic function 

x(Ai,A 2 ,A 3 ) = Ti[\T )(T \D 1 (X 1 )®D 2 (X 2 )®D 3 (X 3 )] 



(0|E/Jl>i(Ai) ® L> 2 (A 2 ) ® D 3 (X 3 )U t \0) 
exp [~(\K\ 2 + \K\ 2 + \K\ 2 ) 
,t 



(6) 



where Xj are complex numbers, Dj(Xj) = exp(Aja] — XjCij) is a displacement operator for 
the j-th mode, and the primed quantities are obtained by using the Heisenberg evolution of 
the modes in Eq.s ((HJ). In formulas 

A'i = fiXi - giX 2 - hiX 3 



A' 2 = -/2A1 + g 2 X 2 + h 2 X 3 



A 3 = -/3Ai+^A 2 + /i 3 A 3 . (7) 



Following Ref.—, the characteristic function can be rewritten as 

1 



x(Ai,A 2 ,A 3 ) = exp 



--x T Cx 
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(8) 



where x T = (xi,X2,X3,pi,p2,P3), (• • -) T denotes transposition, Xj = 2 l ' 2 {p 3 - <.'/<■ 



j = 1, 2, 3, and C denotes the covariance matrix of the Gaussian state, whose explicit 
expression can be easily reconstructed from Eq.s (J7J). The covariance matrix determines the 
entanglement properties of |To). In fact, since |To) is Gaussian the positivity of the partial 
transpose is a necessary and sufficient condition for separability 21 , which, in turn, is deter- 
mined by the positivity of the matrices Tj = AjCAj — iJ where Ai = Diag(l, 1, 1, —1, 1, 1), 
A2 = Diag(l, 1, 1, 1, —1, 1), A3 = Diag(l, 1, 1, 1, 1, —1) and J is the symplectic block matrix 




I being the 3x3 identity matrix. A numerical evaluation of the eigenvalues of Tj shows that 
they are nonpositive matrices Vj. Correspondingly, the state in Eq. (jSJ) is fully inseparable 
i. e. not separable for any grouping of the modes. Notice that the success of a true tripartite 
quantum protocol, as the telecloning scheme described in the following sections, is a sufficient 
criterion for the full inseparability of the state iTo}^ 2 -. 

3. Telecloning of coherent states 

Here we show how the three-mode entangled state described in the previous section can be 
used to achieve telecloning^ of coherent states^, that is to produce two clones at distance 
of a given input radiation mode prepared in a coherent state. Depending on the values of 
the coupling constants of the Hamiltonian (0) the two clones can either be equal one to 
each other or be different. In other words, the scheme is suitable to realize both symmetric 
and asymmetric cloning machines^. This option can eventually be useful to fit the purpose 
of the clones production in order to distribute the quantum information contained in the 
input stat o 29 i 30 i 31 i 32 . Our scheme, which is analogous to that of Ref.— in the absence of an 
amplification process for the signal, is applied to the telecloning of coherent states, whereas 
the state we use to support the teleportation is the three- mode entangled state of Eq. (jHJ). 
For the symmetric case we obtain an optimal cloning machine, achieving the maximum value 
of fidelity allowed in a continuous variable cloning process (F = 2/3) 30 i 31 i 32 . In the case of 
the asymmetric cloning a range of coupling parameters can be found that allows the fidelity 



of one clone to be greater than 2/3, maintaining the fidelity of the other greater than 1/2, 
i.e. the maximum value reachable in a classical communication scheme. 




Fig. 1. Schematic diagram of the telecloning scheme. After the preparation of the state |To) by 
bilinear interactions in a nonlinear crystal (NLC), a conditional measurement is made on the mode 
cti, which corresponds to the joint measurement of the sum- and difference-quadratures on two 
modes: mode a\ itself and another reference mode b, which is excited in a coherent state a, to be 
teleported and cloned. The result of the measurement is classically sent to the parties who want 
to prepare approximate clones, where suitable displacement operations (see text) are performed. 

A schematic diagram of our scheme is depicted in Fig. ^ After the preparation of the 
state | To) a conditional measurement is made on the mode cq, which corresponds to the joint 
measurement of the sum- and difference-quadratures on two modes: mode cq and another 
reference mode b, whose state is to be teleported and cloned. The measurement can be 
described as the following cr-dependent POVM acting on the mode cq 

Ii{(3) = -D{(3)a T D\(3), (9) 

7T 

where D(P), (5 6 C, is the displacement operator, and a is the preparation of b, i.e. the 
state to be teleported and cloned. 

The probability distribution of the outcomes is given by 



P(P) = Tr 123 [|To)(To|n(/3)®I 2 ®I 3 ] 

= -A -V gfffL- + f (p + q \D((3)a T D\(3)\p + q) . (10) 

ir(l + Ni) ^ (1 + Nx)p+Q p\q\ ^ y| K ' ' KH)W Hl V ; 

\ / pq \ j 

The conditional state of the mode a 2 and a 3 after the outcome (3 is given by 
Qp = ^TrjTo><Toin(/?)®I 2 ®I 3 ] 
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P(/J) 7T(1 + JVx) ^(l + iVi)(^+ fc+ ')/2 
x (A: + /| J D(/?V T J Dt(/5)b + g)b,g)(A; ! /|. 




(p + g )!(fc + Q! 
p! g! fc! Z! 



(11) 



After the measurement the conditional state may be transformed by a further unitary op- 
eration, depending on the outcome of the measurement. In our case, this is a two-mode 
product displacement Up — D T ((3) ® D T (f3) where the amplitude (3 is equal to the results 
of the measurement. This is a local transformation which generalizes to two modes the 
procedure already used in the original CV 1 — > 1 teleportation protocol. The overall state 
of the two modes is obtained by averaging over the possible outcomes 



where Tp = Up Qp Up. 

If b is excited in a coherent state a = \z)(z\ the probability of the outcomes is given by 



Moreover, since the POVM is pure also the conditional state is pure. We have Qp = 
with 



i. e. the product of two independent coherent states. The amplitudes are given by 

5 2 p = (z + P)k 2 5 3 p = (z + (3)k 3 , 
where the quantities Kj, j = 2, 3 are given by 





(12) 



|#)) = 1*2/9)2 ® ^3/3)3 , 



(13) 




(14) 



Correspondingly, we have Tp = Up \^p)){{^p\ Ul with 



Up \ijp)) = \zk 2 + P(k 2 - 1)) <8> \zk 3 + /3(k 3 - 1)) . 



(15) 



The partial traces g 2 = Tr 3 [^ 2 3] and g 3 = Tr 2 [f?23] read as follows 



Qj = [ d 2 (3 P Z (P) \zKj + 0( Kj - l)){zKj + P( Kj - 1)1 . 



(16) 
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We see from the teleported states in Eq. (fTtjj) that it is possible to engineer a symmetric 
cloning protocol if N 2 = N 3 = N, otherwise we have an asymmetric cloning machine. 
Consider first the symmetric case. According to Eq.s (j3J) the condition N 2 = N 3 = N holds 
when 

cosnt= 171,2 (17) 

2|72f - iTil 

from which it follows that 



Since |/3 )| 2 = exp{ — 1/3 — /3 | 2 }, the fidelity of the clones is given by 

= \ . (19) 

2 + 3^-2^(2^ + 1) 

As we expect from a proper cloning machine, the fidelity is independent of the amplitude 
of the initial signal and for < iV < 4 it is larger than the classical limit F = 1/2. Notice 
that the transformation Up performed after the conditional measurement, is the only one 
assuring that the output fidelity is independent of the amplitude of the initial state. In Fig. 
12] the behavior of the fidelity versus the average photon number N is shown in the relevant 
regime. We can see that the fidelity reaches its maximum F = 2/3 for N = 1/2 which 
means, according to Eq. ([18)1 . that the physical system allows an optimal cloning when its 



coupling constants are chosen so that I71/72I = v6 — V32 ~ 0.586 . Let us now consider 
the asymmetric case. For N2 7^ A^3 the fidelities Fj = (z\gj\z) of the two clones (|TK|) are 
given by 

F 2 = l - (20) 

2 + iV 3 + 2N 2 - 2^/N 2 (N 2 + N 3 + l) 

F3 = \ = . (21) 

2 + N 2 + 2N 3 - 2^/N 3 (N 2 + N 3 + l) 

A question arises whether it is possible to tune the coupling constants so as to obtain a 
fidelity larger than the bound F = 2/3 for one of the clones, say g 2 , while accepting a 
decreased fidelity for the other clone. Indeed, for example, if we impose F 3 = 1/2, i.e. the 
minimum value to assure the genuine quantum nature of the telecloning protocol, then we 
should choose A^ 3 = \N 2 . In this case the maximum value for F 2 is given by F 2max = 4/5, 
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which occurs for N 2 = 1. More generally, by substituting Eq. (}2*T]) in Eq. pUj) and 
maximizing F 2 with respect to N 2 keeping F3 fixed, we obtain that for N 2 = I/F3 — 1 and 
N 3 = 1/4(1/F 3 - l)- 1 we have 

F 2 = 4 (1 - F3) 
2 (4-3F 3 )' 

which shows that F 2 is a decreasing function of F% and that 2/3 < F 2 < 4/5 for 1/2 < F3 < 
2/3. Notice that the sum of the two fidelities F 2 + F 3 = 1 + 3/4^^ is not constant, being 
maximum in the symmetric case F 2 = F 3 = 2/3. Notice also that the roles of q 2 and £3 are 
interchangeable in the considerations. 

0.7 H ' ' 1 ' ' 1 ' 



^ 0.6- N. 

0.5-1 . 1 , 1 . 1 . — 

12 3 4 

N 

Fig. 2. Fidelity of symmetric clones versus the average (equal) photon number N of modes 
a 2 and 03. 

4. Telecloning with seeded crystal 

In order to confirm the feasibility of the telecloning scheme presented in the previous section 
we now show that the same protocol can be implemented also when the state that supports 
teleportation is generated by Hamiltonian (0) starting from a coherent state in one of the 
modes, rather than from the vacuum. This may be of interest from the experimental point 
of view, since seeding a crystal with a coherent beam is a useful technique to align the setup, 
and allows the verification of the classical evolution of the interacting fields [see Section |U] . 

The analysis of the scheme is analogue to that of the previous Section, however starting 
from the initial state \a, 0, 0) instead of the vacuum. The explicit expression of the evolved 
state |T a ) is derived in appendix El Notice that the conservation law (J2J) implies that the 
populations for seeded crystal Nj a = (T a \ajaj\T Q ) satisfies the relation Aq a — N 2a — N 3a = 
\a\ 2 . We refer the reader to appendix El for the explicit expressions of Nj a and for their 
connections to the populations Nj for vacuum input. 
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A compact expression for the evolved state is the following 

|T«) = D^aM-t)) ® D 2 (-af 2 (-t)) ® D 3 {-af 3 {-t))\T ) , (22) 

where the fj(t), j = 1,2,3 are given in appendix 1X1 Expression (J22J) can be easily derived 
by using the Heisenberg equation of motion for the field-mode a-i(t) (see Eq.s (J5J)). The 
telecloning process proceed as in the previous Section, with calculations performed using 
the shifted Fock basis \i/j n )i = Di(afi)\n)i, and \ip n )j = Dj(—afj)\n)j, j = 2,3. If the 
reference mode b is excited in a pure coherent state a = \z)(z\, then, as in Sectional the 
conditional state is pure gp = \i>p)){{if>p\ with 

|V/9» = |C2^>2 ® |C3/3>3 , (23) 

i. e. the product of two independent coherent states. The amplitudes are given by 
C2/3 = (z + (3 - a/i)«2 - a/ 2 C3/3 = (z + /? - a/i)K 3 - a/ 3 , 

where the quantities Kj, j = 2, 3 are given by Eq. ()14p. The unitary transformation on a 2 
and a 3 that completes the telecloning is now given by 

Up = D\CP - -W2)® D\(fi - k 3 Wi ~ «/s) • (24) 

In fact, the output conditional state coincides with that of Eq. (|15p. so that the partial 
traces are identical to those given in Eq. (|16j). For N 2 = N 3 = N we obtain symmetric 
clones with the same fidelity as in Section El Moreover conditions (fTTj) and (JTHJ) still hold. 
Notice that also the protocol for asymmetric cloning can be straightforwardly extended to 
the present seeded scheme. 

5. Conditional generation of two-mode entanglement 

Another application of the three-mode entangled state of Eq. (J5J) is the conditional gen- 
eration of a two-mode entangled state of radiation by on-off photodetection on one of the 
modes of state |T ). Indeed, it is possible to produce a robust two- mode entangled state 
that approaches a TWB for unit quantum efficiency i] of the photodetector. In the following 
we evaluate some properties of the conditional state when rj 7^ 1 in order to quantify its 
closeness to an ideal TWB. Notice that, due to the well known properties of TWB, this 
scheme also provides a valid check of the whole apparatus from an experimental viewpoint. 
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Let us consider the situation in which a mode of the state |T ), say the third mode, is re- 
vealed by an on-off photodetector. The probability operator measure (POVM) is two-valued 
{n , III}, n + III = I, with the element associated to the "no photons" result given by 

n Q = I x ® I 2 <g> ^(1 - rj) n \n) 33 {n\ . (25) 

n 

The probability of the outcome is given by 

P = Tr 123 [|To)(T |n ] 

= 1 yp( N 2 \ n / N 3 (l - ri) \ m (n + m)\ 
1 + Nt ^ I 1 + ^i / V 1 + JVi J n\m\ 

m,n ' 

= (l + V N 3 )-\ (26) 
while the conditional output state go = 7^Tr 3 [|T )(T | n ] reads as follows 



n-\-n I 

l + r] N 3 / N2 \— {N 3 {l-r))\ m 1 (n + m)\(n' + m)\. , ... 

On = > — t\ ; — ; \n + m,n)(n + m, i 

1 + iVi ^Xl + Nj \ l + Ni J miM n\n'\ 1 A ' 



Remind that iVi = iV2 + iV^. If 77 = 1 this state reduces to the following TWB 

I — 

When the efficiency of the detector is not unitary a question arises on how to quantify the 
closeness of qq to the ideal state |^o)- From an operational point of view, we can evaluate 
the photon number correlation between the first and second mode, which is defined as 

((nx -n 2 ) 2 ) - ((m) - {n 2 )f 
02 " M + M ' (29) 

and is zero in case of TWB. After straightforward calculations we arrive at 

^3(1-^(1 + ^3) 
U2 (l + V N 3 )[2N 2 + N 3 (l - V )\ ' 1 J 

which, for any given value of the quantum efficiency 77, is a decreasing function of N 2 and 
an increasing function of N 3 . A global quantity to characterize the state in Eq. (}2Tj) is the 
fidelity with a reference TWB state. The natural choice for the reference is the TWB \i/)q), 
according to the following argument. At first we calculate the fidelity between state (|27jl 
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and a generic TWB of parameter £ i.e. |£) = -^/l — £ 2 £ n |ra, n.), we have 



.2N 1 + ^3 



1 + 

n + n 



l + iV ly / V 1 + Ni 7 m! V n\n'\ 

i +v n 3 i~e 



2 • 



(31) 



1 + Nl (1-^^/(1 + ^))' 
Then, we look for the parameter £ that maximizes the fidelity. Expression (|3*T|) shows that 



the value of £ maximizing the fidelity is, independently on rj, £ = \JN 2 / (1 + iVi). By 
substituting in Eq. (J3T|) we arrive at 

Ffo) = 1 + ^ 3 . (32) 
w; I + A3 K ' 

Therefore, the maximum fidelity is obtained for rj = 1 and the correct reference state is the 
TWB |^o) • I n conclusion, the state generated through conditional on/off photodetection on 
the third mode of |T ) is a robust two-mode entangled state with a fidelity to a TWB given 
by ()32|). Notice that 77 < F{rf) < 1 for any choice of A3. The same analysis is valid for a 
conditional measurement performed on mode a 2 , in which case we obtain an entangled state 
of modes a\ and 03 [in this case the role of N 2 and A3 should be exchanged in Eq.s (|3U|). 
(13 1|) . and (JH2J)]- On the other hand, we notice that a conditional photodetection on mode 
ai does not lead to an entangled state of modes a 2 and (Z3. 

6. The optical scheme 

An experimental implementation of the scheme proposed in this paper can be obtained by 
using a single nonlinear crystal in which the two interactions described by Hamiltonian 
take place simultaneously. The interactions correspond to two phase-matched second-order 
nonlinear processes in which five fields interact and two of them do not evolve (parametric 
approximation). Among the five fields dj involved in the interactions, a 4 and 05 will be the 
non-evolving pump-fields. 

The energy-matching and phase-matching conditions required by the interactions can be 
written as to± — ujx + u 3 , uj 2 — ^3 + W5, k* = ki + k 3 and k 2 = k 3 + k 5 , being k, the wave- 
vectors (in the medium) corresponding to Uj, which make angles dj with the normal to the 
entrance face of the crystal. It is possible to satisfy these phase-matching conditions with a 
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Fig. 3. Interaction scheme. The pump beams a 4 and a 5 are supposed to impinge on 
the crystal face along the normal. The values of the crystal cut angle, and of the 
interaction angles $2, and $3, are calculated to satisfy the phase-matching conditions. 
The wavelengths of the interacting modes are X(ui) = X(uj 3 ) = 1064 nm, A(co> 4 ) = A^) = 
532 nm and A (0^2) = 355 nm. 

number of different choices of frequencies and interaction angles depending on the choice of 
the nonlinear medium. Here we propose an experimental setup based on a (3— BaB 2 4 crystal 
(BBO, cut angle 32 deg, cross section 10 x 10 mm 2 and 4 mm thickness, Fujian Castech 
Crystals Inc., Fuzhou, China) as the nonlinear medium and the harmonics of a Q-switched 
amplified Nd:YAG laser (7 ns pulse duration, Quanta- Ray GCR-3-10, Spectra-Physics Inc., 
Mountain View, CA) as the interacting fields. We choose a compact interaction geometry 
in which two type I non-collinear interactions with the two pump-beams superimposed in a 
single beam with mixed polarization take place (see Fig. EJ). With reference to Fig. El the 
wavelengths of the interacting modes are \{u\) = A (£^3) = 1064 nm, A(u; 4 ) = A (0*5) = 532 
nm and A (a;2) = 355 nm. The interaction angles, calculated by supposing that the two pump 
beams propagate along the normal to the crystal entrance face, result to be d' = 37.74 deg 
, i?i = —-#3 = 10.6 deg and # 2 = 3.5 deg, and since the crystal we used was cut at 32 deg, 
it had to be rotated to allow phase matching. In order to demonstrate the feasibility of the 
scheme in Fig. 01 we adopted the experimental setup depicted in Fig. EJ The fundamental 
and second harmonic outputs of the Nd:YAG laser were sent to a harmonic separator and 
then each beam was collimated to a diameter suitable to illuminate the BBO crystal. The 
polarization of the second harmonic beam emerging from the laser is elliptic, and the two 
polarization components were separated through a thin-film plate polarizer (Pi in Fig. 0]). 
On the ordinarily polarized component a A/2 plate was inserted to modulate the intensity 
of beam as, without affecting the intensity of the other pump, a 4 . The two beams were 
then recombined through a second thin-film plate polarizer (P 2 ) and sent to the BBO. As a 
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first verification of the effectiveness of the interaction described by the Hamiltonian ((TJ, we 
implemented the seeded configuration discussed in Section 0] by injecting the BBO with a 
portion of the fundamental laser output (see Fig. HJ) to realize the initial condition for field 




S I . 



Fig. 4. Experimental setup. The fundamental and second harmonic outputs of the Nd:YAG 
laser are sent to a harmonic separator and then each beam is collimated to a diameter suitable 
to illuminate the BBO crystal. The polarization of the second harmonic beam emerging from 
the laser is elliptic, and the two polarization components are separated by a thin-film plate 
polarizer (Pi). On the ordinarily polarized component a A/2 plate is inserted to modulate 
the intensity of beam 05, without affecting the intensity of the other pump, a 4 . The two 
beams then recombine at a second thin-film plate polarizer (P 2 ) and are sent to the BBO. 
D 5 and D 2 are pyroelectric detectors. 



As a first quantitative check, we measured the energy, E2, of the beam generated at 002 
as a function of the energy, E 5 , of the ordinarily polarized pump-beam for fixed values of 
the energies of the extraordinarily polarized pump-beam, E4, and of the seed-beam, Ei. 
We preliminarily measured Ei by using a pyroelectric detector (ED200, Gentec Electro- 
Optics Inc., Quebec, QC, Canada) which also allows checking the stability of the source. By 
averaging over more than 100 pulses we found a value of about 48 mJ per pulse, only 50% 
of which is ordinarily polarized, and thus suitable for the interaction. To measure energy 
E 4 we inserted another pyroelectric detector (mod. ED500, Gentec) after P 2 . By averaging 
again over more than 100 pulses we found a value of about 158 mJ per pulse. To obtain 
a reliable measurement of E5 we inserted, on the path of beam as, a cube beam splitter 
and a calibrated glass plate to extract a fraction of the beam. Energy E5 was varied by 
rotating the A/2 plate, and its measurement was performed with the same detector ED500 
as before (see D 5 in Fig.HJ)). To measure the energy E 2 of the output pulses we used another 
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pyroelectic detector (PE10, Ophir Optronics Ltd., Jerusalem, Israel, see D 2 in Fig.HJ)). The 
values of E5 and E2 were measured simultaneously as averages over the same 20 laser shots 
at each rotation of the A/2 plate. In Fig we show the measured values of E 2 (open circles), 
as a function of the measured values of E 5 . We can compare the experimental results with 
the field evolution calculated according to the classical equations 24 



u 2 CiE 4 ■ c 2 E 5 



cos 



(yc 2 E b - dE 4 z^j - 1 



E 1 



(33) 



wi (c 2 E 5 - C1E4) 

where c\ = 8.3 x 10 4 (Jm 2 ) -1 and c 2 = 2.6 x 10 5 (Jm 2 ) -1 are the coupling constants that apply 
to the present interactions. In Fig we show the values (full circles) of E2 as calculated 
according to Eq. (}3*3*j) for the experimental values of E 5 , and for fixed values Ei = 24 mJ 
and E4 = 158 mJ. The agreement between measured and calculated values is excellent. 




Fig. 5. Comparison of the experimental results with the field evolution calculated according 
to equation Open circles: measured values of the energy of field a 2 as a function of 

the measured values of the energy of pump-field a 5 . Full circles: values of the energy of 
field a 2 as calculated from the classical evolution of the interacting fields as a function of 
the measured values of the energy of pump- field a^. 



7. Conclusions and outlooks 

We have suggested a scheme to generate fully inseparable three-mode entangled states of 
radiation based on interlinked bilinear interactions taking place in a single nonlinear 
crystal. We have shown how the resulting three-mode entanglement can be used to realize 
symmetric and antisymmetric telecloning machines that achieve optimal fidelity for coherent 
states. An experimental implementation involving a BBO nonlinear crystal is suggested and 
the feasibility of the scheme is analyzed. Preliminary experimental results are presented: 

15 



as a first quantitative check, we measured the energy of the beam generated at U2 as a 
function of the energy of the ordinarily polarized pump, for fixed values of the energies of 
the extraordinarily polarized pump-beam, and of the seed-beam. The agreement between 
measured and calculated values is excellent. 

To realize the telecloning protocol described in Section |H] we need to generate three-mode 
entanglement from vacuum. This should be possible by implementing the same experimental 
setup as in Fig. 01 with a different laser source able to deliver a higher intensity. In fact, 
we plan to use a mode- locked amplified Nd:YLF laser (IC-500, HIGH Q Laser Production, 
Hohenems, Austria) with which it is easy to achieve an intensity value of 50 GW/cm 2 in a 
collimated beam. Since such a value was enough to generate bright twin beams in a 4-mm 
thick BBO crystal, it should allow us to obtain the three-mode entangled state described in 
this paper, not only by seeding the crystal but also for initial vacuum state. 
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Appendix A: Heisenberg evolution of modes 

In this section we calculate the dynamics generated by the Hamiltonian (|T} in the Heisenberg 
picture. The equations of motion are given by 



a 2 = -i^ 2 a 3 



(Al) 



This system od differential equations can be Laplace transformed in the following algebraic 
system 

a\(0) + fJ>a\(fi) = Z7 1 5 3 (/i) 
a 2 (0) + //a 2 (/i) = -z'720 3 (/i) 
03(0) + /xa 3 (/x) = -z'7ial(/i) - ry 2 a 2 (/x) , 



where we have defined the Laplace transform of a^t) 

a,(/x) = / dt e~^aj(t) . 
Jo 

The determinant of the system (|A2|) is 

A = M o + r)( M -r), 



where T = a/|7i| 2 — | ^2 1 2 ? therefore its solution reads 

40") = ^ (l72| 2 + /i 2 )ol(0) +7i7 2 °2(0) +«7 1 /ia 3 (0) 
a 2 (/i) = ~t~ f-7i72Oi(0) + (/i 2 - |7i| 2 )a 2 (0) - ij 2 ^a 3 (0) 



03(1") = ^ -z'7i/ia}(0) -i/i7 2 a 2 (0) +// a 3 (0) 



The solution of system (jAl|) follows from anti-transforming Eq. ()A5|) . We have 

a\(t) = fia\(0) + f 2 a 2 (0) + f 3 a 3 (0) 
a 2 (t) = gia\(0) + g 2 a 2 (0) + g 3 a 3 (0) 
a 3 (t) = h ia \(0) + h 2 a 2 (0) + h 3 a 3 (0) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 
(A7) 
(A8) 
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where the coefficients are given by 



fi(t) = ^[|7ii 2 cosm-| 72 | 2 ] (A9) 

= [cos fit -1] (A10) 

/ 3 (t) = sin (fit) (All) 

9i(t) = ^[1- cos fit] (A12) 

= ^O7i| 2 -I7 2 | 2 cosfit] (A13) 

= sin (fit) (A14) 

^(t) = S i n (fit) (A15) 

/ia(t) = _i2i s in(fit) (A16) 

/i 3 (t) = cos (fit) (A17) 



and fi = iT = yj 



72 - 7i 



2 



Appendix B: Schrodinger evolution in a seeded crystal 

In this appendix we derive the explicit expression of the evolved state from \a, 0, 0). We can 
write the Hamiltonian (jlj as follows: 

H int = 7l K^+^J + h.c, (Bl) 

with the definitions K = a\a^ and J = 02(13. To calculate the evolved state we can proceed 
by factorizing the temporal evolution operator of the system; to this purpose we introduce 
the following operators 

Ji = a\a\ + 0303 , J2 = a\a3 — c4°2 , M = a\a^ , 

which form with K and J a closed algebra. Actually, the temporal evolution operator can 
be written in the following way: 

fjM _ e (3iK^ e faM^ e foJ^ e PAJi e foJ2 e PeJ e (hK e f3 6 M ^ ( B2 ) 

which allows us to calculate the evolution of a generic initial state as a function of (3+. In 
the case under investigation we obtain: 
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E~ or 
-=|n,0,0> 
vn! 

M,o> 

v^(ae ft ) n , 



e -^-e /3l ^ t e /32Mt e /33jt e /34Jl V s 



-'n! Vp' Vn! 



n.p 



in \/( n +P + 9) ! i 



It can be demonstrated- that 



e"^ 4 £ ^(ae^r v ^7^^> + p + g,p, q) . (B3) 



v^TivT ' ^ V i + JVi ' V i + JVi ' 

Moreover, for the population with initial vacuum Nj = (T |ajaj|T ) and initial seed Nj a 
(T a \ajaj\T a ) we have the relations 

* = ^z!#, iv 2 = ^, iv 3 



1 + M 2 ' l + |a| 2 ' J 1 + lal 2 ' 
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